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PART  ONE 


THEORY  AND  EXAMPLES 


T .  INTRODUCTION 

Formulas  for  the  computation  of  plane  wave  transmission  through 
a  rectangular  aperture  in  a  perfectly  conducting  plane  are  derived  in 
Part  One.  The  computer  programs  which  use  these  formulas  are  given  in 
Part  Two.  The  general  theory  of  solution  is  derived  in  a  previous  re¬ 
port  r 1 1 .  Basical lv,  the  procedure  is  an  application  of  the  method  of 
moments  to  an  integral  equation  formulation  of  the  problem.  The  unknown 
to  he  determined  is  the  equivalent  magnetic  current  over  the  aperture 
rtgion,  which  is  proportional  to  the  tangential  electric  field  in  the 
aperture.  The  solution  is  expressed  in  terms  of  an  aperture  admittance 
matrix,  which  is  dual  to  the  impedance  matrix  for  the  complementary  con¬ 
ducting  plate.  Once  the  equivalent  magnetic  current  is  obtained,  the 
electromagnetic  field  can  be  computed  via  potential  integrals.  The  nota¬ 
tion  used  in  this  report  is  the  same  as  that  used  in  [1].  We  abstract 
equations  from  this  previous  work  ar.  we  need  them,  referring  to  them  by 
equation  number.  We  do  not  attempt  to  summarize  the  theory  here. 

i  •  1 ous  studies  of  aperture  problems  include  those  for  small  aper¬ 
tures  am'  rinse  for  circular  apertures  f4,5].  Some  results  for 


[1]  R.  F.  Harrington  and  J.  R.  Mautz,  "A  Generalized  Network  Formulation 
for  Aperture  Problems,"  Scientific  Report  No.  8  on  Contract  F19628-73-C- 
0047  with  A.F.  Cambridge  Research  Laboratories,  Report  AFCRL-TR-75-0589 , 
November  1975. 

[2]  H.  A.  Bethe,  "Theory  of  Diffraction  by  Small  Holes,"  Phys.  Rev., 
vol.  66,  pp.  163-182,  October  1944. 

[3]  C.  J.  Bouwkarap,  "Diffraction  Theory,"  Repts.  Progr.  in  Phys., 
vol.  17,  pp.  35-100,  1954. 

[4]  G.  Bekefi,  "Diffraction  of  Electromagnetic  Waves  by  an  Aperture  in 
a  Large  Screen,"  .Toum.  Appl,  Phys.,  vol.  24,  No.  9,  pp.  1123-1130, 
September  1953. 

[5]  C.  Bouwkamp,  "Theoretical  and  Numerical  Treatment  of  Diffraction 
Through  a  Circular  Aperture,"  IFF.F  Trans,  on  Antennas  and  Propagation, 
vol.  AP-18,  No.  2,  pp.  152-176,  March  1970. 
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aperture  in  a  conducting  plane 


apertures  of  arbitrary  shape  have  been  obtained  using  Babinet's 
principle  plus  a  wire  grid  approximation  to  the  complementary  con¬ 
ducting  plate  [6,7].  The  reader  may  consult  these  papers  for  other  j  > 

references. 

IT.  STATEMENT  OF  THE  PROBLEM 

Figure  1  shows  the  prohlem  to  be  considered  and  defines  the  coordi¬ 
nates  and  parameters  to  be  used.  The  infinitely  conducting  plate  covers 
the  entire  z=0  plane  except  for  the  aperture,  which  is  rectangular  in 
shape  with  side  lengths  L^Ax  and  L^Ay  in  the  x  and  y  directions,  respec¬ 
tively.  The  excitation  of  the  aperture  is  a  uniform  plane  wave  incident 
from  the  region  z  <  0.  The  field  to  be  computed  is  the  far  zone  magnetic 
field  in  the  region  z  >  0,  at  the  angles  0,<j>. 

The  solution  is  expressed  in  terms  of  the  equivalent  magnetic  current 
M  =  z  x  E,  where  z  is  the  unit  z-directed  vector  and  E  is  the  electric 
field  in  the  aperture.  To  compute  M,  we  use  a  linear  expansion  in  terms 
of  basis  functions  M.  and  evaluate  the  coefficients  by  the  method  of 
moments.  This  involves  determining  a  generalized  admittance  matrix, 
evaluated  in  Section  III,  and  an  excitation  vector.  To  determine  the 
field  produced  by  M,  we  need  a  measurement  vector.  The  excitation  and 
measurement  vectors  for  the  present  problem  are  of  the  same  form,  and  are 
evaluated  in  Section  IV. 

III.  ADMITTANCE  MATRIX 

According  to  [1,  Eq.  (28)]  and  [1,  Eq.  (10)],  the  admittance  matrix 
fY]  is  given  by 

Vi  =  (y3  +  Yb)i]  "  '  4<V  (1) 

[6]  A.  T.  Adams,  C.  B.  Vamado,  D.  E.  Warren,  "Aperture  Coupling  by 
Matrix  Methods,"  1973  IEEE  EMC  Symposium  Record,  New  York  City, 

June  1973,  pp.  226-240. 

f 7 ]  J-L  Lin,  W.  L.  Curtis,  M.  C.  Vincent,  "On  the  Field  Distribution 

of  an  Aperture,"  IEEE  Trans,  jn  Antennas  and  Propagation,  vol.  AP-22, 

No.  3,  pp.  467-471,  May  1974. 
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where 


V  •  WJ 

<~1 


We  choose  the  set  of  testing  functions  W  equal  to  the  set  of 
expansion  functions  .  The  rectangular  aperture  0  x  Lx^x» 

0  <_  y  <_  LyAy  where  and  Ly  are  integers  is  divided  into  rectangular 

subareas  of  length  Ax  in  x  and  Ay  in  y.  The  set  of  expansion  func- 
tions  is  split  into  a  set  of  x  directed  magnetic  currents  and  a  set 

y 

of  y  directed  magnetic  currents  defined  by 

p  =  1,2 , . . ,L  -1 

mx  lur  r,  -  *  tx(x)  py(y)  »  x  < 

~P+(q-l)(Lx-l)  -p  q  *  q  =  1,2, ...L 

y 


~P+(q-l)hx  ^  l  T^(y)  PX<x)  , 


P  =  1 *2 *  * • • 
q  =  l,2,...Ly-l 


where  TX(x)  and  T^(y)  are  triangle  functions  defined  by 
P  q 


TX(x)  = 
P 


x  -  (p-l)Ax 
Ax 

(p+l)Ax  -  x 
Ax 


(p-1)  Ax  <_  x  _<  pAx 


pAx  <_  x  <_  (p+l)Ax 


|  x  -  pAx  I  >_  Ax 


(q-l)Ay  <_  y  <_  qAy 


qAy  <_  y  (q+l)Ay 


|y  -  qAy|  >_  Ay 


and  PX(x)  and  P^(y)  are  pulse  functions  defined  by 
P  q 
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Px(x)  = 


P^(y)  = 


1  (p-1 ) Ax  x  <  pAx 


0  all  other  x 


1  (q-l)Ay  y  <  q Ay 


0  all  other  y 


XV  x  y 

The  magnetic  charge  sheets,  sav  p  and  p7  associated  with  M  and  M 
are  obtained  from  (5)  as 


(Px(x)  -  P*  . (x) )Py (y) 

p  p+i  q 


p+ ( q - 1 ) (Lx-1) 


p+(q-l)I- 


(Py(v)  -  Py  , (y))Px(x) 

_ 9 _ a±i _ e _ 

-JwAy 


X  V 

Introduction  of  the  two  types  of  expansion  functions  IT  and  M7  and 
the  two  types  of  testing  functions  and  M7  into  (8)  gives  rise  to 
four  Y  submatrices  defined  by 


Yil  "  4*“ 


(Fj  •  Mj  +  4y“)ds 


V  V 

where  u  is  either  x  or  y  and  v  is  either  x  or  y.  In  (18),  and 
are  the  electric  vector  and  magnetic  scalar  potentials  due  to  . 

The  integrations  over  the  "field"  magnetic  current  and  charge 

and  p^  explicit  in  (18)  are  approximated  by  sampling  the  integrands 
1 

at  two  points.  Hence, 


Y™  “  A^'AXAy  [J  (^1  '  *}x  ,y  +  I  (l]  *  ..  v 

J  p  q  P+1.' 


1  ,„V 


ImAx  (^)xp,yq  +  1o!ax  (^)xp+1,yq1 
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VU  •  '2  'll  •  i>Vyq  *  *  ->  ’  -  VVi 


1  „  +  77X7  -  1 


ju»Ay  .1  xp,yq  >Ay  j  *p’Vi 
where  v  is  either  x  or  y  and 


(20) 


x  =  (p  -  . 5) Ax 
P 

y  =  (q  -  .5)  Ay 

q 

To  determine  p  and  q  in  terms  of  i  in  (19)  or  (20),  refer  to  (10) 
or  (11). 

Substitution  of  (10),  (11),  (16),  (17),  (3)  and  (A)  into 
(19)  and  (20)  yields 


(21) 

(22) 


yXx  =  jAxAX  jl  x  (s  t-q)  -  J  1  (s-p+1,  t-q) 
ij  nn  2  c  ^  x 


+  (s~P-  ±-2121  I  (s-p+l,  t-q)  +  \  I  (s-p-1,  t-q)  -  tezBz2Z2)  T  (s-p-1, t-q) 
2  c  z  x  ^ 

+  -~-r  (1  (s-p+l,  t-q)  -  21  (s-p,  t-q)  +  r  (s-p-1,  t-q))]  (23) 

k  Ax 


/V*  =  1- 

unk2 


I  (s-p,  t-q)  +  I  (s-p+l,  t-q) 

C  C 


+  T  (s-p,  t-q-1)  -  1  (s-p+l,  t-q-1) ] 

c  c 


(24) 


% 


rnk 


I  (s-p,  t-q)  +  I  (s-p,  t-q+1) 

c  c 


+  I  (s-p-1,  t-q)  -  I  (s-p-1,  t-q+1)] 

C  C 


(25) 
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vIj  -J^r 'iv-”’  “<■>-  '-*«> 


+  —  ^3/2)  tc(s-p,  t-q+^)  +  j  Iy(s-p,  t-q-1)  - 


(t-q-3/2)  .  , 

- i -  ICs-p.t-q 

l  c 


+  — — -X  (T  (s-p,  t-q+1)  -  21  (s-p,  t-q)  +  I  (s-p,  t-q-1))]  (26) 

,  L  .  L  C  c  c 

k  Ay 

where  n  ~  ~  376.730  ohms  Is  the  Intrinsic  Impedance  for  empty  space 

and  where 


(t+l/2)Ay  <s+l/2)Ax  _jk  J~2^2 

|  dy  [  dx  — — ~  (27) 

J  J  r 2  2 

v=(t-l/2)Av  x=(s-l/2)Ax  Vx  +  y 


Ic(s,t)  =  k 


(t+l/2)Ay  (s+l/2)Ax 


T  (s,t)  =  ~~ 
x  Ax 


-jk  Jx2+ y2 


J  1  2  2 

y=(t-l/2)Ay  x=(s-l/2)Ax  vx  +  y 


(t+1/2) Ay  (s+1/2) Ax  _jR 

ly{s't)  =  a7  _  ydy  j  dx  e  ■■  . 

y=(t-l/2)Ay  x»(s-l/2)Ax  vx  + 


I  (s,t)  is  even  in  both  s  and  t ,  I  (s,t)  is  odd  in  s  and  even  in  t, 
c  x 

and  I  (s.t)  is  even  in  s  and  odd  in  t.  In  (23)  to  (26),  is  the 

y  v  ^ 

interaction  of  the  expansion  function  M  with  the  testing  function 

u  x  J 

M  .  For  the  testing  function  M  , 

**si 


i  =  P  +  (q-l)(Lx-l), 


p  =  1,2 , . . .L  -I 
x 


q  =  1 , 2  . . .  L 


whereas  for  the  testing  function  M^, 


1  =  P  +  (q-l)I-  , 


f  P  “  1,2,. ..Lx 
x*  l 

i  q  =  1,2, ...L-1 
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1  -  L  <  (s-p)  <  I,  -2 

x  —  —  X 

2  -  L  <  (t-q)  <  L  -1 

y  —  —  y 

but  because  (2 A)  is  odd  about  (s-p)  «  -1/2  and  odd  about  (t-q)  =  1/2, 


(s-p)  =  0,1 ,2,. . .L  -2 
(t-q)  =  1,2,3, ...L  -1 


(35) 


is  sufficient,  Tn  (25) 


1 


t 


2  -  I.  <  (s-p)  <  L  -1 

X  “  X 

f 

1  -  L  <  (t-q)  <  L  -2  f 

y  —  —  y 

hut  because  (25)  is  odd  about  (s-p)  *  1/2  and  odd  about  (t-q)  =  -1/2, 


(s-p)  =  1,2,3,. ..Lx-1 
(t— q)  =  0,1,2, ...L  -2 
is  sufficient.  Finally,  in  (26), 

1  -  L  <  (s-p)  <  L  -1 

x  —  —  x 

2  -  L  <  (t-q)  <  L  -2 

y  -  y 

but  because  (26)  is  even  in  both  (s-p)  and  (t-q). 


(s— p)  =  0,1,2,... 1 
(t-q)  =  0,1,2,... L^— 2 

is  sufficient.  From  inspection  of  (23)  to  (26)  and  (34)  to  (37), 


(37) 


s  —  —  1,0,1,... 1 


t  -  —  1 , 0 , 1 , . . . L^— 1 


(38) 


is  adequate  in  (27)  to  (29) . 

The  integrals  (27)  to  (29)  are  evaluated  by  using  the  following 
four  term  approximation 


-jkr  .2 

e"j  r»  e  °[1  -  jk(r-rQ)  -  —  (r-v  ■  6  •-  0- 


-r  )2  +  (r-r  )3]  (39) 


where 
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r 


(40) 


rQ=  ~J (sAx)2  +  (tAy)2  . 
Substitution  of  (39)  into  (27)  yields 


Ic(s,t)  =  [k(l  +  jkro 


,2  2  J  3 
k  r  Ik  r 
o  o. 


dxd> 

r 


„  2  2 

?  jk  r 

+  k  (-j  +  krQ  +  — j 


Ml 


,  i  Jkr 

dxdy  +  k  (-  - - ~) 


(41) 


rdxdy 


+ 


2 

r  dxdy]e 


(42) 


where  the  limits  on  all  the  integrals  in  (42)  are  the  same  as  those 
in  (27).  The  approximations  to  (28)  or  (29)  are  given  by  (42)  with 
an  additional  factor  of  either  —  or  in  the  integrands.  Three  of 
the  required  integrals  are 


dxdy  =  AxAy 

(43) 

xdxdy  =  sAx2Ay 

(44) 

2 

ydxdy  =  tAxAy 

(45) 

The  indefinite  integrals  associated  with  the  rest  of  the  required 
integrals  are 


II 

II 


x  log(y  +  r)  +  y  log  (x  +  r) 


rdxdy  -  log(y  +  r)  +  ^-  log 


(x  +  r) 


(46) 


(47) 
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r  dxdy  -  -X^  — 

(48) 

OfL-V-  +^log  (y  +  r) 

(49) 

2  2  4 

xrdxdy  =  yr(yj  +  ^-)  +  log  (y  +  r) 

(50) 

,  -  2  2 

xr  dxdy  =  x  y(^—  +  ^-) 

4  0 

(51) 

-  Y  +  V  lo*  (x  ♦  r) 

(52) 

2  2  4 

yrdxdy  «  xr(jy  +  lo8  (x  +  r) 

(53) 

,  ,22 
yr  dxdy  =  y  x  (^-  +  ^-) 

(54) 

The  reader  can  verify  (46)  to  (54)  by  showing  that,  in  each  case,  the 


mixed  second  partial  derivative 


3x3y 


of  the  right  hand  side  is  equal 


to  the  integrand  on  the  left  hand  side.  The  definite  integral  is  ob¬ 
tained  from  the  indefinite  integral  by  adding  the  indefinite  integral 
evaluated  at  both  upper  limits  to  that  at  both  lower  limits  and  sub¬ 
tracting  both  evaluations  of  the  indefinite  integral  at  the  mixed  (one 
upper,  one  lower)  limit3. 


IV.  PLANT  WAVE  EXCITATION  AND  MEASUREMENT  VECTORS 

The  plane  wave  excitation  vector  of  [1,  Eq.  (32)]  and  the 
plane  wave  measurement  vector  of  [1,  Eq.  (37)]  are  of  the  same  form 
except  for  a  minus  sign.  We  therefore  need  to  evaluate  only  one  of 
them,  say  the  measurement  vector  We  specialize  it  to  four  princi¬ 

pal  plane  patterns  as 


1 
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(P™)0y  =  _  2  [j  Mj  .  0  ejkx  C0S  6  dxdy  (55) 

■pert. 

,„mu.  .  [I  wu  -  Ikx  cos  6  .  . 

(P.  )  «  -  2  M  •  y  eJ  dxdy  (56) 

i  yy  J  J  '*'1  -* 

apert. 

<pr>4,  -  -  *  \\  «?  ■  i  *lky  “■  *  “v  <”> 

apert. 

,_mu.  .  ff  ,.u  -  jky  cos  ,  ,  /co. 

(P^  )  =  -  2  ill  *  £  e  dxdy  (58) 

apert . 

The  superscript  u  is  necessary  because  has  been  split  up  into 

and  of  (10)  and  (11).  In  (55)  to  (58),  9,  y,  <f>,  and  x  are  unit 

vectors  in  the  0,  y,  and  x  directions  respectively  where,  as  shown 

in  Fig.  1,  0  is  measured  from  the  positive  x  axis  in  the  y  *=  0  plane 

and  $  is  measured  from  the  positive  y  axis  in  the  x  =  0  plane.  For 

measurement  vectors,  0°  0  <_  180°,  O'1  <_  $  ".80°.  (P™U)g  is  for  a 

0  polarized  measurement  in  the  y  =  0  plane,  (P^)^  is  for  a  y  polarized 

measurement  in  the  y  =  0  plane,  (P^U)^x  is  for  a  $  polarized  measurement 

in  the  x  =  0  plane,  and  (p™U)xx  is  for  a  x  polarized  measurement  in  the 

x  =  0  plane.  Because  our  set  of  testing  functions  W  is  the  same  as  the 

+i 

set  of  expansion  functions  M  ,  the  plane  wave  excitation  vector  P  of 
[1,  Eq.  (32)]  is  obtained  by  putting  180°  <  6  <  360°,  180°  ^  360°  in 

the  negative  of  one  of  the  equations  (55)  to  (58). 

Substituting  (10)  and  (11)  into  (55)  to  (58)  we  obtain,  with  the 
help  of  [9] 

[9]  H.  B.  Dwight,  Tables  of  Integrals  and  Other  Mathematical  Data, 
fourth  edition,  Macmillan  Co.,  New  York,  1961,  Eq.  567.1. 
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l 
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V.  REPRESENTATIVE  COMPUTATIONS 

A  versatile  computer  program  has  been  developed  using  the  pre¬ 
ceding  formulas.  This  program  is  described  and  listed  in  Part  Two 
of  this  report.  Some  representative  computations  obtained  with  this 
program  are  given  in  this  section. 


The  first  computations  were  made  for  a  narrow  slot,  of  width 
A/20  and  of  variable  length  L.  The  far-zone  quantity  plotted  was  the 
transmission  cross  section,  defined  as  [1,  Eq.  (39)] 

t  -  2iTr^|H  (67) 

m' 

where  H  is  the  component  of  magnetic  field  being  considered.  We 
m 

use  the  notation: 

9  9 

x.  =  2rr  |naj  in  the  y  =  0  plane, 

0y  9  (68) 

2  2 

t  =  2irr  Ih  I  in  the  x  =  0  plane, 
xx  1  x 1 

For  the  case  being  considered,  the  orthogonal  components  of  H  in 
these  two  planes  were  zero.  Figure  2  shows  plots  of  T0y  and  txx  for 
x-directed  slots  of  width  A/20  and  length  (a)  L  =  A/4,  (b)  L  =  A/2, 

(c)  L  =  3A/4,  and  (d)  L  =  A.  In  all  cases  the  excitation  was  due  to 
a  plane  wave  normally  incident  on  the  conducting  plane  with  the  mag¬ 
netic  field  in  the  x  direction.  Note  the  large  transmission  cross 
seo^ion  for  L  =  A/2,  case  (b),  due  to  the  slot  being  near  resonance. 

The  plots  of  t  are  of  the  same  form  as  scattering  cross  section  from 
the  complementary  conducting  strips,  as  known  from  Babinet's  principle. 

Figure  3  shows  plots  of  the  equivalent  magnetic  current  in  the 
aperture  region  for  the  same  slots.  Since  M  =  £  *  £,  they  are  also 
plots  of  the  tangential  component  of  E  in  the  slots.  Again  note  the 
large  value  of  M  for  the  case  L  =  A/2,  which  is  near  resonance.  Note 
also  that,  for  short  slots  (L  <_  3A/4),  the  M_ is  almost  equiphasal  and 
closely  approximated  by  a  half  sine  wave. 

Next,  computations  were  made  to  test  the  rate  of  convergence  of 
the  solution  as  the  number  of  subsections  was  increased.  A  slot  of 
width  A/10  and  length  2A  was  chosen  for  the  study.  Again  the  excitation 
is  a  plane-wave  normally  incident  on  the  conducting  plane  with  the  mag¬ 
netic  field  in  the  x  direction.  Figure  4  shows  plots  of  and 
for  the  cases  (a)  39,  (b)  19,  (c)  9,  and  (d)  4  triangular  expansion 
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Transmission  cross  section  for  slots  of  length  L  in  the  x 
direction  and  width  A/20  in  the  y  direction.  (a)  L  =  A/' 
(b)  L  ■  A/2,  (c)  L  *  3A/A,  (d)  1  *  A.  Excitation  is 
by  a  plane  wave  normally  incident  on  the  conducting  plane 
with  magnetic  field  in  the  x  direction. 


MAGNITUDE  MAGNITUDE 


ansmission  cross  section  when  the  number  of  expansion  furc 
ons  is  (a)  39,  <b)  19,  (c)  9,  and  (d)  4.  r°mP^ations  are 
r  a  slot  of  length  2X  in  the  x  direction  and  width  X/10  in 
e  y  direction.  Excitation  is  by  a  plane  wave  normally 
cident  on  the  conducting  plane  with  magnetic  field  in  the 

A  4  f*  •(  nn  . 


magnitude  magnitude 


PHASE 


Transmission  cross  sections  for  a  square  aperture  of  side 
length  L,  excited  by  a  plane  wave  with  H  in  the  xz  plane 
and  incident  at  an  angle  0  from  the  normal  direction  in 
the  H  plane.  (a)  L  =  X/4,  9=0.  (b)  L  =  X/4,  0  =  45°. 


functions  respectively.  Note  that  the  patterns  (a)  and  (b)  are  essen¬ 
tially  the  same,  and  pattern  (c)  is  only  slightly  different.  Tliev  dif¬ 
fer  appreciably  from  (d) ,  which  results  from  only  4  expansion  functions. 

The  difference  in  the  solutions  as  the  number  of  expansion  functions  is 
decreased  is  better  illustrated  hv  plots  of  M,  as  shown  in  Fiji.  5.  These 
are  for  the  same  cases  as  the  corresponding  cases  of  Fig.  4.  It  can  be 
seen  clearly  how  the  computed  equivalent  current  in  the  slot  region  changes 
as  the  number  of  subsections  is  reduced.  As  a  rule  of  thumb,  for  near¬ 
field  quantities  (such  as  M)  one  should  use  subareas  of  length  >/10  or  less 
and  for  far-field  quantities  (such  as  t)  length  J/5  or  less. 

Finally,  Fig.  6  shows  some  computations  for  wider  apertures  and 
excitations  by  waves  not  normally  incident  on  the  conducting  plane.  All 
cases  shown  are  for  square  apertures,  of  side  length  L.  Figures  6(a)  and 

(b)  are  for  1,  =  >/4,  with  the  plane  wave  normally  incident  for  (a)  and 
incident  45°  from  the  normal  direction  in  the  H  plane  for  (b) .  Figures 
6(c)  and  (d)  are  for  I,  =  >/ 2,  with  the  plane  wave  normally  incident  for 

(c)  and  45°  from  the  normal  direction  in  the  H  plane  for  (d).  Note  that, 
for  the  relatively  small  slots  chosen,  there  is  little  difference  in  the 
shapes  of  the  patterns  as  the  incident  wave  direction  is  changed  from  the 
normal  direction.  There  is,  however,  an  appreciable  difference  in  the 
amplitudes  of  the  patterns. 

VI.  DISCUSSION 

The  computer  program.  Part  Two,  is  written  explicitly  for 
rectangular  apertures,  but  the  formulas  are  valid  for  any  aperture 
composed  of  rectangular  subsections.  Other  apertures,  such  as  I. -shaped, 
T-shaped,  square  O-shaped,  etc.,  could  be  treated  by  appropriately 
changing  the  computer  program.  Apertures  of  arbitrary  shape  could  he 
treated  hv  approximating  them  by  rectangular  subsections.  As  with  all 
moment  solutions,  the  size  of  the  apertures  which  can  be  treated  depends 
upon  the  size  of  the  matrix  which  can  be  computed  and  inverted.  The  ex¬ 
amples  Indicate  that  the  rectangular  subsections  should  have  side  lengths 
not  greater  than  0.2  wavelengths  for  reasonable  accuracy. 
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The  aperture  admittance  matrix  has  application  to  any  problem 
in  which  one  region  is  hounded  by  a  plane  conductor,  as  shown  in 
reference  [1],  Hence,  it  can  be  used  for  waveguide-fed  apertures  in  a 
ground  plane,  and  for  cavity-backed  apertures  in  a  ground  plane.  It  is 
planned  to  treat  these  latter  two  problems  in  future  reports. 
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PART  TWO 


COMPUTER  PROGRAMS 

l.  DESCRIPTION  OF  THE  MAIN  PROGRAM 

The  main  program  computes  the  complex  coefficients  V  which 

n 

determine  the  magnetic  current  ^according  to  [1,  Eq.  (5)],  the  trans¬ 
mission  coefficient  T  of  [1,  Eq.  (44)],  and  four  patterns  of  the  trans- 

2 

mission  cross  section  [1,  Eq .  (40)]  per  square  wavelength  t/X  .  The 
2 

four  patterns  of  ~[/\  are  written  on  the  first  record  of  direct  access 
data  set  6.  The  main  program  calls  the  subroutines  LINEO,  YMAT,  and 
PLANE  which  are  listed  later  on  in  this  report. 

One  data  card  is  read  early  in  the  main  program  according  to 

READ  (1,11)  LX,  LY,  LI,  NTH,  OX,  DY,  TH 
11  FORMAT  (413,  3E14.7) 

The  L^  and  1,  appearing  in  (10)  are  read  in  through  LX  and  LY  respec¬ 
tively.  Here,  LX  >_  2  and  LY  ^  1  which  means  that  the  long  dimension  of 
a  rectangular  aperture  only  one  subsection  wide  must  lie  along  the  x 
axis.  DX  is  Ax/X  and  DY  is  Ay/X  where  Ax  and  Ay  appear  in  (12)  and  (13) 
and  X  is  the  wavelength.  The  plane  wave  excitation  vector  [1,  Eq.  (32)] 
is  the  negative  of  expression  (54  +  LI)  where  LI  is  either  1,2,3  or  4  and 
where  the  angle  (either  0  or  <p)  in  degrees  appearing  in  equation  (54  +  LI) 

is  TH.  The  four  patterns  of  the  transmission  cross  section  [1,  Eq.  (40)] 

2 

per  square  wavelength  t/>  are  generated  by  evaluating  the  plane  wave 
measurement  vectors  (55)  to  (58)  at  angles  (0  or  if)  equal  to  (J-l)*180./ 
(NTH-1)  degrees,  .1  =  1,2, ...NTH. 

Minimum  allocations  are  given  by 

COMPLEX  Y (N*N) ,  P(4*N),  V(N) 

DIMENSION  TA(4*NTH) 

where 

N  =  (I.X-1)  *LY  +  LX*  (LY-1)  (69) 
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A 


Statement  27  uses  LX,  LY,  DX,  and  DY  to  store  — — . —  [Y3  +  Y^ ] 

jAxAy 

3  b 

where  [Y  +  Y  ]  is  the  admittance  matrix  appearing  in  [1,  Eq.  (14)] 

by  columns  in  Y.  Here,  n  =  =  376.730  ohms  is  the  intrinsic  impedance 

for  empty  space.  At  the  time  statement  27  is  executed,  LX  and  LY  are 
still  the  original  input  data,  namely  the  numbers  of  subdivisions  in  x 
and  y  but  [IX  and  DY  are  kAx  and  kAy  instead  of  the  original  input  data 

Ax/'  and  Ay/X. 

Statement  28  inverts  the  N  by  N  matrix  stored  in  Y. 

Statement  2b  uses  TH,  LX,  LY,  DX,  and  DY  to  store  - — —  P^  in 

2AxAy  m 

P(r  +  (V-1)*N),  K  =  1,2, 3, 4,  m  =  1,2,...N,  where  P1  is  the  plane  wave 

m 

excitation  fl,  Fq .  (32)]  and,  in  particular,  the  negative  of  expression 
(S4  +  K) ,  At  the  time  statement  29  is  executed,  TH  is  the  angle  0  or  4> 
(see  (55)  to  (58)  and  Fig.  1)  in  radians  which  specifies  the  direction 
from  which  the  exciting  plane  wave  comes,  LX  and  I,Y  are  the  numbers  of 
subdivisions  in  x  and  y,  DX  is  kAx  and  DY  is  kAy. 


Nested  DO  loops  16  and  17  multiply  the  matrix  stored  in  Y  by  the 
column  vector  stored  in  P(1  +  (I,I-1)*N)  through  P(LI*N)  and  use  the  con¬ 
stant  UV  =*  -j 2 -m ri  in  order  to  store  V  of  [1,  F.q.  (14)]  in  V.  Substituting 
fl,  F.q.  (28)]  into  [1,  Eq.  (13)]  we  obtain 

[Yhs]V  =  j  T1  (70) 

which  simplifies  [1,  Fq.  (44)]  to 


T  = 


2qS  cos  0 


.  -  ->•!  * 

Re  (V  I  ) 


(71) 


inc. 


-  M  * 

-VI 

DO  loop  16  accumulates  ; —  in  U2 .  Statement  31  stores  the  trans- 

K  2 AxAy 

mission  coefficient  T  of  (69)  in  T. 


[1,  Fq.  (40)]  simplifies  to 

,  4 


t/X2  = 


3  2 
32  m  n 


IP^I 


.ti2 


(72) 


1 


24 


DO  loop  19  stores  the  transmission  cross  section  per  square  wavelength 
2 

t/A  of  (72)  in  TAU(K).  Statement  30  uses  TH  =  ( J-l) *u/ (NTH-1)  radians, 

LX,  LY,  DX,  and  DY  to  store  — ?m  in  P(n  +  (K-1)*N).  For  the  ~  P1” 

2AxAy  n  2AxAy 

piny 

stored  in  P(1+(K-1)*N)  through  P(K*N),  DO  loop  21  accumulates  — — - —  in  Ul. 

2AxAy 

Next,  i /A 2  of  (72)  is  stored  in  both  TAU(K)  and  TA(J+(K-1) *NTH) . 

Statement  32  writes  TA  on  the  first  record  of  data  set  6  for  possible 
input  to  the  plot  program  listed  later  on  in  this  report. 
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C.  LISTING  OF  THE  MAIN  PROGRAM  AMD  SAMPLE  OAT  A 

r 

//  EXEC  WATFIV 

//GP.F T06F001  00  DSNAMF  =EE0034.RFVl  , 0  I S P  =  OL D , UN  I T= 33 30 , 

//  Of  B= (RFC  FM  =  VS ,  RLKS I  Z F= 2 596 , L P FCL  =  25 92  ) 

/ /GO .SYS  IN  00  * 

iJOR  MAUTZfTI ME=1*PAGES=40 

MAIN  PROGRAM 

f  THIS  PROGRAM  CALLS  THE  SUBROUTINES  L I NE 0 , Y MA T , PL ANF 

COMPLEX  U,UV, YI2500) ,P(200) ,U1,U2,V(50) ,CONJG 
0  I  ME  NS  I  ON  T Ab!4) , T A<  1168) 

PI =3. 141503 
ET A=376. 730 
U* (0.,  1.  ) 

UV=2.*PI *ETA*U 

R  E  AO  I l *  1  1 1  LX ,LY, L l ,NTH ,DX  ,DY ,TH 

11  FORMAT (413, 3F14. 7  ) 

WP I TF<  3,  12  J  L  X  ,LY ,Ll  ,NTH, OX  ,nY, TH 

12  FORMAT  I  •  LX  LY  L  I  N TH • , 5 X  ,  • DX •  , 1 2 X , • D Y • , 1 2X , • TH • /I X , 41 3 , 3E 1 4 . 7) 
RK  =  2  •  *  P  l 

OX  =0  X*  RK 
OY=DY*RK 
P8=180./PI 
TH  =  TH/ P8 

27  CALL  YMAT (LX,LY,OX,DY,Y) 

WP ITF ( 3,  13)  (  Y<  I ) , I =1 ,3  I 

1 3  FORMAT  (  •  Y* /<  1X,6E  11.4)  ) 

N-  (LX-1)  *LVHX*(L  Y-l  ) 

28  CALL  L INEOt  N, Y) 

WRITE  (3, 13)  (Yd),  1=1,3) 

20  CALL  PLANE! TH, LX, LY.DX ,DY,P) 

WR  ITE  (3,  14M  PI  I),  1  =  1,3) 

14  FORMAT ( •  P' /IX, 6E  1  1.4) 

!A=l+(LI-l)*N 

I B  = I A  +  N-l 
U2  =0. 

DO  16  J=  1 , N 
U1  =0. 

J1  =  J 

DO  17  I  =  I  A  ,  IP 
Ul=Ul*V( Jl) *P( I ) 

J  1  =  J  1  ♦  N 

17  CONTINUE 
V( J ) =U1*UV 
J 1 = I A+J-l 

U2=U2«-V(J)*CCNJG!  P(J1)  ) 

16  CONTINUE 

WR I T  E ( 3 , 24 )  (V(  I),  1  =  1, N) 

24  FORMAT!'  COEFFICIENTS  V  OF  MAGNETIC  CURRENT  EXPANSION  FUNCTIONS' 
1/1 1X.6EL1.4) ) 

31  T  =  REAL(U2)/«LX*LY*ETA*SIN!TH)  ) 

WR  I  TE ( 3  ,  18)  T 

18  FORMAT!'  TRANSMISSION  COEFFICIENT  T=',E14.7) 

CT=DX*DY/< P I ♦ E  T  A) 

CT=CT*CT/(8.*PI ) 

DT  H  =  P  I  / ( NTH- 1  ) 

WP  I  T E ( 3 , 23 ) 

23  FORMAT  l ' 0  ANGL E • , 4 X, • T AU  l ' , 7X , • TAU2 • , 7X , • T AU 3 • , 7X , * T A U4'  ) 

DO  19  J= 1 »  N  TH 


TH  = ( J-  1 ) *DT  H 

30  CALL  PLANE! TH,LX, LY,DX , OY ,P I 
TH  =  T  H*  P8 
J1=0 
J2  =  J 

00  20  K=  1  *  4 
U  1  =  0  . 

DO  2 1  1=  1,N 
J  1  =  J  1  ♦  1 

UUUUPI  J1)*V(  I  ) 

21  CONTINUE 

H=  U 1 *C  CN JG ( U  1  I 

TAU(K)=CT*H 
TA( J2»=TAU( KJ 
J2=J2+NTH 
20  CONT I NUE 

WR I T  E ( 3 , 22 )  TH, (T AU( I)  ,  1= l ,4) 

22  FORMAT! IX, F 7.2, 4E  1  l.4| 

19  CONT INUE 

KA=J2-NTH 
REWlNO  6 

32  WRITE  (6)  (TA(J»  ,J=1,KA> 

STOP 

FND 

t  DA  T  A 

8  1  119  0.5000300F -01  0 . 5030030E -0 1  0. 2 700030E+03 

t  ST  HP 
/* 

// 

PRINTED  OUTPUT 

LX  LY  LI  NTH  OX  DY  TH 

5  1  1  19  0. 5  OOOOOOF -  01  0. 5000000F -0 1  0. 2 700000E+03 

Y 

-0.  153  1F+02-0.  6525F-01  0 . 664  bF  +  0  1  -0 . 6  46  3  E -0  1  0.  1  31  2  F  +  01  -0 . 6  2  7 5  E -0  1 

Y 

-0.9999F-01  0.4015F-02- J.6822F-01  0 .553 7F-02-0 . 5 1 3 1 E-01  0.5523E-02 

p 

-0. 1000E+01-0. 30  3  4F - 06- 0. IODOE+01-0.6067F-06-0.1000F+01-0 .9 10  IE -06 
COEFFICIENTS  V  OF  MAGNETIC  CURRFNT  FXPANSION  FUNCTIONS 
0.451  IF  4-02  0.5916F+03  0.6238E+02  0.8153E+03  0.6238E+02  3.8153E+03 
0.451  1  F  >02  0.5916F+03 
TRANSMISSION  COEFFICIENT  T=  0  .  I  l  4 1 2  54  E  +  00 

ANGLE  TAU1  TAU2  TAIJ3  TAU4 

0.00  0.  OOOOF  *00  0 • 00  OOF  +  00  O.OOOOF+OO  0.2186F-02 
10.00  0. 5885F-04  O.OOOOF+OO  O.OOOOF+OO  0.2186F-02 
20.00  0.2308F-03  C.OOOOE+OO  O.OOOOF+OO  0.21B8E-02 
30.00  0.  50 16E-0  3  O.OOOOE+OO  O.OOOOF+OO  0.2I90F-02 
40.00  0. 8462F-0  3  O.OOOOE+OO  O.OOOOF+OO  0.2193E-02 
50.00  0. 1 228F-02  O.OOOOF+OO  O.OOOOF+OO  0.2196F-02 
60.00  0. 1602F-02  O.OOOOF+OO  O.OOOOE+OO  0.2199E-02 
70.00  0.  1 9 1 8E-02  O.OOOOF+OO  O.OOOOF+OO  0.  2202F-02 
80.00  0.2129E-02  O.OOOOE+OO  O.OOOOF+OO  0.2203F-02 
90.00  0.2204F-02  C.OOOOF+OO  O.OOOOF+OO  0.2204E-02 
100.00  0.2129E-02  O.OOOOE+OO  O.OOOOE+OO  0.2203F-02 
110. 00  0. 1918E-02  O.OOOOF+OO  O.OOOOE+OO  0.2202E-02 
120.00  0. 1602F-02  O.OOOOF+OO  O.OOOOF+OO  0.2199F-02 
130.00  0. 1228F-02  O.OOOOF+OO  O.OOOOE+OO  0.2196F-02 
140.00  0 .  84 62 E -0 3  O.OOOOE  +  OO  O.OOOOF+OO  0.2193F-02 

27 


150.00  0. 5016E-03  O.OOOOE  +  OO 
160.00  0. 23  08E-0  3  O.OOOOE+OO 
170.00  0.5885E-04  O.OOOOE  +  OO 
180.00  0. 7663E- 1 5  0.0000E+00 


O.OOOOE+OO  0.2190E-02 
0. OOOOE  +00  0. 2 188E-02 
0. OOOOE+OO  0.2186E-02 
0.  OOOOE -*-00  0.2186F-02 


C  LISTING  OF  THE  SUBROUTINE  LINEO 

f 

SUBROUTINE  LINEOILL.CI 
COMPLFX  CI2500) ,STOR,STO,ST,S 
DIMENSION  L  R ( 50) 
on  20  1  =  1,  LL 

lp  m  =  i 

20  CONT  INUE 
Ml  =  0 

DO  18  M=l,Ll 
K  =  M 

K2  =  M  UK 

SUAPSIPEAL  (C(K2)  )  UABS(AIMAG(C(K2)  )) 
DO  2  I =M, LL 
K 1  =  M  U  l 

S2  =  APS(REAL  (C(K1)  )  UABS  (  A  I  MAG  (C  (K  1 )  )) 
I F ( S  2- SI )  2,2,6 

6  K=  I 

S 1  =  S2 

2  CONTINUE 

L  S  =  L  R  C  m  ) 

L  P  (  M  )  =  l  R  (  K  ) 

LR ( K  )  =LS 
K2  =  MUK 
STOP  =C ( K2) 

J1  =0 

00  7  J=l ,LL 
K1  =  JUK 
K2  =  J  1  +M 
ST  0=  C ( K  1 ) 

C ( K  1  )  =C ( K2 ) 

C( K2 ) =STO/STCR 
JUJHU 

7  CONTINUE 
K 1  =  M  U  M 

CIK1  )  =  1 . /ST  OR 
DO  ll  I  =  1 , L  L 
IF(I-M)  12,11,12 
12  K  1  =  M 1  ♦  I 
ST=C(Kl) 

C  (  K  1 )  =  0. 

J1  =  0 

Dn  10  J=1,LL 
K  1  =  J  U  I 
K2  =  J 1  4-M 

C( K1 ) =C(K1»-C(K2) *ST 
Jl  =  JULL 
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10  CONTINUE 

11  CONTINUE 
Ml  =M1HI 

1  8  CONTINUE 

J1  =0 

DO  9  J=1,LL 
IM J-LRI J»  )  14,8,  14 

14  l R  J  =  L  R ( J  1 

J?=(LPJ-l)*LL 

2  1  DO  10  1  =  1 , l  L 

K2  =  J?M 
K  1  =  J  1  ♦  I 
S  =  C  <K2) 

C ( K2 ) *C ( K  1  > 

C(  K 1  I =S 
13  CONTINUE 

LR ( J) =LR(L8 J) 

L  R  <  L  P  J  »  = I  RJ 

IF ( J-IRI J» )  14,8,  14 

8  J1=JULL 

9  CONTINUE 
RETURN 
END 
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II.  DESCRIPTION  OF  THE  SUBROUTINE  YMAT 


The  subroutine  YMAT(I.X,  LY,  DX,  DY,  Y)  uses  the  input  variables 
LX,  LY,  OX,  and  DY  to  calculate  and  to  store  by  columns  in  Y  where 

[Y]  is  the  admittance  matrix  (1)  dealt  with  in  Part  One.  With  regard 
to  the  input  variables,  IK  and  LY  are  the  numbers  Lv  and  L  (see  (10))  of 
subdivisions  in  x  and  y,  DX  is  kAx  and  DY  is  kfiy  where  Ax  and  Ay  (see  (12) 
and  (13))  are  the  lengths  of  the  subdivisions  in  x  and  y.  We  require  that 
LX  2.  and  l-'1’  >_  1. 

Minimum  allocations  are  given  by 

COMPLEX  TC (.11)  ,  TX(J1),  TY(J1),  YXX(J2),  Y (N*N) 

.11  =  (LX  +  1)  *  (LY  +  l) 

.12  =  MAX  ((IK-1)  *  LY.  LX  *  (LY-1)) 

N  =  (TJX  -1)  *  LY  +  LX  *  (LY-1) 

Here,  MAX  denotes  the  maximum  value. 

Nested  DO  loops  15  and  16  put  I^(s,t)  of  (42)  in  TC(JST),  l^(s,t) 
in  TX(.rST),  and  r^(s,t)  in  TY  (JST)  where 

s  =  .IS  -  1 
t  =  JT  -  1 

JST  =  s  +  2  +  (t+l)(Lx  +  1) 

As  mentioned  in  Part  One,  the  expressions  for  I  (s,t)  and  1  (s,t)  are 

x  y 

similar  to  (42).  The  variables  x  and  y  of  integration  in  (43)  to  (54) 
are  changed  to  kx  and  ky  with  the  result  that  x  and  y  is  replaced  by  kx 
and  ky  everywhere  on  the  right  hand  sides  of  (43)  to  (54)  and  the  dangling 
factors  of  k,  k2,  k3,  and  k4  in  (42)  disappear.  The  logic  inside  DO  loop  16 
is  best  understood  by  building  up  a  table  of  variables  in  YMAT  versus  expres¬ 
sions  in  *-erms  of  variables  appearing  in  Part  One,  Section  III. 


30 


Variables  In  YMAT 


Expressions  in  Part  One,  Section  III 


f 


YL 

(t  -  . 5) kAy 

! 

Yl! 

(t  +  .5) kAy 

1, 

| 

XL 

(s  -  .5)kAx 

1 

XU 

(s  +  .5)kAx 

1 

R1 

kro  22  33 

k  r  k  r 

U1 

i  ,  -i  o  0 

1  +  jkro  '  2  "  J  6 

.2  2 
k  r 

U2 

-j  +  krQ  +  j 

kr 

U3 

1  .  o 

~  2  J  2 
-)kr 

o 

e 

EX 

SI 

kx  log  (y  +  r)  +  ky  log(x+r) 

evaluated  at  x,y  limits 

S3 

k^x^  log(y+r)  +  k^y^  log(x+r) 

,3  „  ,3  3  ,33 

evaluated  at  x,y  limits 

S5 

-k-P  +  Y  log  (v+r)  + 

log(x+r)  evaluated  at 

x,y  limits 

TC  (JST) 

I  (s,t) 
c 

.  3  2  4 

S5 

k  +  +  log(y+r)) 

,  9  2  2 

evaluated  at  x,y  limits 

S6 

V  x  v  (“T"  +  evaluated  at 

4  6 

x,y  limits 

TX  (JST) 

Ix(s,t) 

,  3  2  4 

S5 

k  (y§ —  +  +  log(x+r)) 

_  _  2  2 

evaluated  at  x,y  limits 

S6 

5  2  y  x 

k  y  x  (^7—  +  evaluated  at 

4  o 

x,y  limits 

TY  (JST) 

ly(s,t) 

i 
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In  the  preceding  table,  the  first,  second,  and  third  S5  and  the  first 
and  second  86  correspond  respectively  to  the  first,  second,  and  third 
uses  of  85  and  the  first  and  second  uses  of  S6  in  DO  loop  16. 

If  I  =  1,  the  logic  between  the  statements  44  and  18  uses  the 

fact  that  (27)  is  even  in  s  to  store  T  (-1,0)  in  TC(L  +2).  Similarly, 

c  x 

I  (-1,0)  is  stored  in  TX(L  +2)  and  I  (-1,0)  in  TY(1.  +2).  If  L  4  1,  the 
x  x  y  x  y 

logic  between  statements  44  and  18  stores  T  (s,-l)  in  TC(s+2), 

c 

s  =  -  1,0,1  ... .I_x-1  and  Tc (-1 , t)  in  TC(1  +  (t+l)(I.  +  1)),  t  =  0,1, ...L  -1 

and  similarly  for  I  and  T  . 

x  y 

Nested  00  loops  19  and  20  store  -rr~ — YXX  where  YXX  is  given  by  (23) 

j AxAy  ij  ij 


with 


(s  -  p)  =  .IS  -  2 
(t  -  q)  =  -IT  -  2 


in  YXX(s-p+l  + (t-q) (I.  —1) ) .  See  (34)  for  bounds  on  (s-p)  and  (t-o). 

Nested  00  loops  24,  23,  22,  and  21  store  — — YXX  where  i  and  j 


are  given  by  (30)  and  (32)  where 


jAxAy  ij 


p  =  JP 
q  =  JQ 
s  =  JS 
t  =  JT 

in  Y(i  +  (j-l)*N)  where  N  is  given  by  (69).  If 

s  -  p  0 
t  -  q  >_  0 

the  subscript  for  YXX  Inside  nested  00  loops  24,  23,  22,  and  21  is 

s  -  p  +  1  +  (t-q)  (I.  -  1) 

x 

The  more  general  subscript 


|s  -  pi  +  1  +  |  t-q  |  (I-x  -  1) 


is  a  consequence  of  the  fact  that  Y^j  of  (23)  is  even  in  both  (s-p)  and 
(t-q) . 
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Nested  no  loops  25  and  26  store 
(24)  with 


yyx 
jAxAy  ij 


where  Y'* 


Is  given  by 


s  -  p  -  JS  -  2 
t  -  q  =  JT  -  2 

in  YXX(s-p  +  1  +  (t-q-l)(Lx  :  1)).  See  (35)  for  bounds  on  (s-p)  and  (t-q) . 

Nested  DO  loops  30,  29,  28,  and  27  store  —  yT*  where  i  and  1  are 
f  ’  jAxAy  ij  J 

given  by  (31)  and  (32)  where 


p  =  JP 
q  =  JQ 
s  =  JS 
t  =  JT 

in  Y ( (L  -  1)  L  +  i  +  (j-l)*N).  If 

x  y 

s  -  p  >_  0 

t  -  q  >_  1 

the  subscript  for  YXX  inside  nested  DO  loops  30,  29,  28,  and  27  is 

s  -  p  +  1  +  (t-q-1)  (L^  -  1) 

The  more  general  subscript 

I*  -  P  +  jl  +  J+  (Jt  -  q  -  \\  ~  |)ax  -  1) 

for  YXX  is  a  consequence  of  the  fact  that  Y^*  of  (24)  is  odd  about 

1  1 
s-p  =  -  —  and  odd  about  t-q  = 

tt  n  xv  xv 

Nested  DO  loops  31  and  32  store  — — - —  Y.,  where  Y,,  is  given  by 

r  jAxAy  ij  ij  n 

(25)  with 

s  -  p  -  JS  -  2 

t  -  q  -  JT  -  2 

in  YXX  (s-p  +  (t-q)(L^  -  1)).  See  (36)  for  bounds  on  (s-p)  and  (t-q). 

Nested  DO  loops  36,  35,  34,  and  33  store  Y^  where  i  and  j 

are  given  by  (30)  and  (33)  where 
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in  Y (N* (I,  -  1)  *  L  +  i  +  (j-l)*N).  If 

x  y 

s  -  p  >  1 
t  -  q  >_  0 

the  subscript  for  YXX  in  nested  DO  loops  36,  35,  34,  and  33  is 

s  -  p  +  (t-q)  (L^  -  1). 

The  more  general  subscript 

!  s  "  P  -  jl  +  J+(|t-q  +  |-|  -  j)  (hx  ~  1) 

for  YXX  is  a  consequence  of  the  fact  that  Y*T  of  (25)  is  odd  about 

1  1 
s-p  =  —  and  odd  about  t-q  =  -  y. 

TTfi  yy  yy 

Nested  DO  loops  37  and  38  store  •  Y*  .  where  Y"  is  given  by 

.1  uxfiy  ij  ij 

(26)  with 

s  -  p  =  JS  -  2 
t  -  q  =  JT  -  2 

in  YXX  (s-p  +  1  +  (t-q)  L  ) .  See  (37)  for  bounds  on  (s-p)  and  (t-q). 

777-1  yy 

Nested  DO  loops  42,  41,  40,  and  39  store  —  Y-  where  i  and  \ 
are  given  by  (31)  and  (33)  where 

p  =  JP 
q  =  -TQ 
s  =  JS 
t  =  JT 


in  Y ( (N+l)  *  (L  -  1)  *  I.  +  1  +  (1-1)  *  N).  Tf 


s  -  p  _>  0 
t  -  q  >_  0 

the  subscript  for  YXX  in  nested  DO  loops  42,  41,  40,  and  39  is 
s  -  p  +  1  +  (t-q)L 
The  more  general  subscript 

!s-p[+l+|t-q[Lx 

for  YXX  is  a  consequence  of  the  fact  that  Y^  of  (26)  is  even  in 
both  (s-p)  and  (t-q). 
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C  LISTING  OF  THE  SUBROUTINE  YMAT 

C 

SUBROUTINE  YM AT  I L X ,L Y , OX , OY  ,  Y  ) 

COMPLEX  U,Ul,U2,U3 ,U4, EX, TCI  100),  TX1 100) , TYI  100) ,YXX(  100) , Y12500) 
0X2= 1 . /(OX*OX| 

OY?=l./(OY*OY) 

DXDY=DX*DY 
NX=ILX-1)*LY 
NY= I LY-1 ) *L  X 
N=NX+NY 
L  X  P  =  L  X  f  1 
LYP=L Y+l 
LXM=LX-l 
LYM=L Y-l 
U= (0. , l. ) 

U4  =  .  I  66666 7 *U 

JST=LX«-l 

DO  15  JT=1,LY 

JST=JST+ l 

YL=< JT-1.51 *  0  Y 

YU=YL+OY 

YL2= YL*YL 

YU2= YU*YU 

Y1  =  I JT-1 ) *D Y 

Y2  =  Y  l*Yl 

on  16  J  S  =  1  ,  L  X 

XL= ( JS-l .5 ) *0X 

XU  =  X  L  *-D  X 

XL2=XL*XL 

XU2=XU*XU 

X  1  =  ( JS-l  )*DX 

X2  =  X  1  *  X  1 

R2=  X2*Y2 

R l =S0PT(P2) 

Rl)  1=  1 .  5* R  2 

U1  =  RIJ1  +R  1*  (  1  1666667*R2  )  *U 

U2  =  R l-RU  1*U 
U3=- . 5-. 5*R l*U 

EX=cns(Rll - U  *  S I N I Rl) 

JST= JST> 1 
R5=XL2+YL2 
R6=  XU2+YL2 
R  7  =  X  L  2  +YU2 
R8  =  XU2+YU2 
P1=SQRT(RS) 

R2=SQRT( R6) 

R3=SQPT( R7) 

P  4  =  S  OR  T  (  R  8  ) 

AYL=YL*AL3G(  (  XU+R  2  )  /  (  X  L  +R  1  )  ) 

AYU=YU*ALOG(  (  XU  +  P  4  J  /  (  XL  «-R  3  I  ) 

AXL  =  XL*ALOG  (  I  YU*P  3)  /  IYL«-R1  )  ) 

A  X  U=  XU*A  LOG  I  <YU  +  P4)/(YL«-R?I) 

S1=AXU-AXL+ AYU-AYL 
AYL= YL*AYL 
AYU= YU*A YU 
AXL=XL*AXL 
A  XU  =  XU* A  XU 

S3  =  XU* AXU-XL*AXL*  YU*AYU-YL* AYL 
XY 1  =  XL  *YL 
XY2  =  XU*YL 
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XY3-XL*YU 
XY4=  XU*YU 

S5  =  .  333333  3*  (  X  Y4*  R  4-  XY  3*  R  3  -XY  ?*  R?  <•  X  Y  1  *R  1  )  ♦  .  I  66666  7*  S3 
rr  (  JS  r  )  =  <S  l  *Ul«-OX  DY*U2«-S5*U3  +  .3  33  3  3  33*(  XY4*R  8-XY3*R  7-XY2*R6*X  Y  1*R5 
l  ) *  U4 ) *  E  X 
VR  l  =  YL*Rl 
YR2=YL*P2 
YR  3  =  YU*R  3 
YR  4  =  YU*R4 

S6=.8333333E-l*<YR4*R8-YR3*R7-YR2*R6+YRl*P5)+.125*(XU2* I YR4-YR2) -X 
L  L  2* ( YR3-YR 1  ) *XU2* A  XU- XL  2* A  XL  ) 

S6=.25*0Y*<  XU2* XU  2 -XL  2* XL  2 )♦.  33  333  33*  X1*DX*(  YU2* YU- YL 2* YL ) 

TX(JST)=.5*  (  YR4-YR  3-YR2+YR  1*  A  XU- A  XL  )*Ul  +  Xl*0X0Y*U2«-S5*U3+S6*J4 

TX ( JST ) =  T  X ( JST )*E X/OX 

XR 1  =  XL*R  l 

XR  2=  XU*R2 

XR3=XL*R3 

XR  4=  XU*P4 

SS=.83333J3E-1*(XR4*RB-XR  3*R7-XR2*P  6+  XR  1  *R  5 )  ♦  .  1 2 5* (  YU 2*  (  X R4 -X R  3>  - Y 
1 L2* ( XR2-XR1  )  *YU2*  AYU-YL2*AYL  1 
S6-.25*0X*( YU2*YU2-YL2*YL2)*«  33  33  333*Y1*DY*(  XU2*XU-XL2*XL) 

TY( JST  )=.5*( XR4-XR  3-XR2+XP  1  ♦ A YU-AYL  ) * U 1 ♦ Y 1 *D XOY* U2* S 5*U 3+ S6* J4 
T  Y<  JST  )  =  TY<  JST) *E  X/DY 
16  CONTINUE 

16  CONTINUE 
IF(lYM)  44,44,45 

44  J1=LXP»1 
J2  =  J  1 *2 

TCI Jl)  -  T  C  I  J  2  I 
TX ( j 1 ) =-TX ( J 2  > 

TY ( Jl ) =TY( J2) 
r,r  TO  46 

45  Jl =2*L  XPtl 

OP  17  JS  =2 , L  X  P 
J  l  =  J  I  ♦  1 

rr i js ) =tc< j  1 1 

TX ( JS) =TX( j  1  ) 

TY( JS)=-TY( Jl) 

1 7  CONT I NUE 
Jl=  1 

nr  i8  jt=i,lyp 

J  2  =  J  1  ♦  2 

Tf (  J  1  )  =  T  C  (  J  2  ) 

TXI  J  1)  =-  TX  (  J  2  ) 

TYI Jl  ) =  T  Y ( J  2  I 
Jl =J 1 *LXP 

1 8  CONT INUE 

46  J4  =  L  X ♦ 2 
J  Y  =  0 

DO  19  JT  =2 , L  Y  P 
00  20  JS  =  2 , L  X 

J  3  =  J4 
J4  =  J  4  ♦  1 
J5  =  J  4+  1 
JY  =  J Y*  1 

YXX(JY)  =  ,5*(TC(J4)»(JS-.5)*TC(J5)-US-3.5)*TC(J3)-TX(J5)*TXU3H*D 
1 X2*(TC IJ5)-2.*TC(  J4)»TC(J31  I 
20  CONT I NUE 
J4  s  J  4 ♦ 2 

19  CONTINUE 


JY  =  0 
DO  24  JT  =  1 1 l Y 
DP  23  JS  =  l »  L  X  M 
DO  22  J0=1,LY 
JTQ  =  LXM*IABS(JT-J0KI 
DO  21  JP=1»LXM 
Jl  =  J  T  Q ♦ I AflS (JS-JP  ) 

JY= JY+l 

YI  JY )  =  YXX<  J  1  ) 

21  CONTINUE 

22  CCNT 1  NOE 
J  Y  =  J  Y  ♦  NY 

23  CONTINUE 

24  CONTINUE 

IE ( LY^.EQ.O )  RETURN 
J4=2*l XP+1 
JY  =  0 

DO  25  JT=3,IYP 
DC  26  J S  =  2  t  L  X 

J  Y  =  J  Y  ♦  1 
J4  =  J4+  1 
J3= J4-LXP 

YXX ( JY )=  <-TC< J4)»TCI J3  »  *TC  « J4  +  1 >-TC( J3*l>) /DXDY 

2  h  r  PNT  I  MJF 
14  =  J  4  ♦  ? 

2  5  r  -SNT  IMlf 

IV  rN  X 

on  30  JT  =  1  * l y 
DP  29  JS=ltLXM 
DO  2  8  JO  =  l,LYM 
JTQ=?*< JT-JQJ-1 
J  2  -  L  X  M*  (  IA«S(JT0»-l)/2 
DO  27  JP=1,IX 
JSP=2* (JS-JP) *1 
J1 =J2* ( I A8S( JSP»+ 1 1/2 
I Y  =  J  Y  ♦  1 

Y<  JY) =YXX<  J l ) 

IF ( JTO+jSP.tT.OI  Y  <  J Y ) =- Y { JY ) 

27  CONTINUE 

28  CONTINUE 
JY= JY+NX 

29  CONTINUE 

30  CONT INUF 
JY  =0 

J4  =  L  XP  +  2 
DO  31  JT  =  2  t  L  Y 
DO  32  JS=3,IXP 
J3  =  J4 
J4  =  J  4  +  1 
J5=J4*IXP 
J  Y  =  J  Y  ♦  1 

YXXI  JY)  =  (-TC(  J4H-TC<  J3I  ♦TCI  J5)-TC  (  J  5- 1 )  >/DXOY 
32  CONTINUE 
J4  =  J  4  +  2 

31  continue 

j Y  =  N*  N  X 

DO  36  JT=1,LYM 
00  35  JS  =  l *  L X 
DO  34  J0= 1 , L Y 
JTQ=2*(  JT-JQH-1 


J2=LXM*(  I A8  S  C JTQI-D/2 
DO  33  JP=1,LXM 

jy=jy*i 

JS P  =  2* (JS-JP)-l 
Jl=J2M!  ABSIJSPl+l  1/2 
Y ( JY )=YXX(J1) 

IF< JTQ*JSP.LT.O)  Y(JY)=-yjY> 

33  CONTINUE 
‘  34  CONTINUE 

JY= JY+NY 
35  CONTINUE 
I  36  CONTINUE 

JY  =  0 
J4=l X+2 
DO  37  JT=2  *  L  Y 
DO  38  JS=2,LXP 
JY= JY*1 
J  4  =  J  4*  I 
J5  =  J4  +  LXP 
J3  = J4-LXP 

YXXI  JY»  =  .5*ITC(J4»*IJT-.5)*TCU5)-(  JT-3.5)*TC(  J3I-TYI  J5H-TYI  J3M  +D 
IY2*ITC(J5)-2.*TC(J4)*TC(J3I* 

38  CONTINUE 
J4= J4+ 1 

37  CONTINUE 
J Y= ( N+ 1 I *N  X 
DO  42  JT=1,LYM 
00  41  JS  =  1 t  LX 
DO  40  JQ=1,LYM 
JTQ  =  LX*IARS(JT-JQ1-H 
DO  39  JP  =  1 1  L X 
Jl  =  JTQ«-I  AB5  (  JS-JP  ) 

JY=JY*  1 
Y(JYI=YXX( Jll 

39  CONTINUE 

40  CONTINUE 
JY= J Y+NX 

41  CONTINUE 

42  CONTINUE 
RETURN 
END 


l 


i 
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III.  DESCRIPTION  OF  THE  SUBROUTINE  PLANE 


The  subroutine  PLANE (TH,  LX,  LY,  DX,  DY,  P)  uses  input  variables 
TH,  LX,  LY,  DX,  and  DY  to  store  times  the  plane  wave  measurements 

(59)  to  (66)  in  P(i),  P(N  +  i) ,  P(N  +  i),  P(N  +  N  +  i),  P(2N  +  i) , 

X  X 

P(2N  +  N  +  i) ,  P(3N  +  i) ,  and  P (3N  +  N^  +  i)  respectively  where 


i  =  p  +  (q-l)(Lx  -  1) 
i  =  p  +  (q-1)  L^ 


N 

X 


(L 


15  Ly 


N  =  N  +  L  (L  -  1) 
x  x  y 


in  (59)  and  (65) 
in  (62)  and  (64) 


Both  angles  A  and  4  appearing  in  (59)  to  (66)  are  equal  to  TH  radians. 

The  arguments  LX  and  LY  of  PLANE  are  the  numbers  L  and  L  of  subdivisions 

x  y 

in  the  x  and  y  directions  and  DX  and  DY  are  the  electrical  lengths  kAx  and 
kAy  of  the  x  and  y  subdivisions.  We  require  that  LX  2,  and  LY  ^  1. 


Nested  DO  loops  81  and  87  store 


1 

2AxAy 


times  (59)  with 


p  =  JP 

q  -  JQ 

in  P(p  +  (q-1)  (L^  -  1)).  Nested  DO  loops  82  and  88  store  times  (62) 

with 

p  =  JP 
q  =  JO 

in  P(N  +  +  p  +  (q-l)I.x).  Nested  DO  loops  83  and  84  store  2&x&y  t^mes  (®4) 

with 

p  =  JP 
q  =  JQ 

in  P(2N  +  N^  +  p  +  (q-Dl.^).  Nested  DO  loops  85  and  86  store  2AxAy  times  (65) 
with 

p  =  JP 

q  -  JQ 

in  P(3N  +  p  +  (q-1) (I.  -  1)). 


i 


! 


I 


t 
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listing  of  thf  surroutinf  planf 

SUP  ROUT  I NF  PI ANFI TH, LX , L  r«  OX, DY  »  P ) 
COMPLEX  U,U1.P<200I 
U= (0.  ,  1. ) 

LXM=L  X-l 
LYM=LY-l 
NX=LXM*L  Y 
N=NX+-LYM*LX 
N4=N*4 

00  89  J=l.N4 
P( Jl =0. 

89  CONTINUE 
SN  =  S I  N  (  T  H) 

CS=COS<TH) 

X2  =DX*CS 

X3=.S*X2 

Sl=  —  SINf X  3  »  t X3 

S2  =  S  I*  S1*SN 

00  81  JP=lfLXM 

S5=JP*X2 

Ul  =  S2*(COS(S5I>U’*,SIN(SS)  ) 

J1=JP 

DO  87  JQ=1,LY 
P( J1I=UI 
Jl =J1+LXM 

87  CONTINUE 

81  CO NT  I NUF 
IFILYM.EO.OI  GO  TO  90 
00  82  JP=1,LX 

S9  = ( J  P- «  9 ) * X  2 

U1 =51* (COS< S5) +U*SIN(S5»  ) 

J 1  =  N*  NX*  Jp 

00  88  JQ=1,LYM 

P ( J  1 ) =U 1 
J1=J1+LX 

88  CONTI NUF 

82  CONTINUE 
90  Y2  =  DY*C  S 

Y3= . 9*Y2 
S1=-SIN( Y3J /Y3 
S2=Sl*Sl*SN 
jl  =  ?*N«-NX 

IFILYM.EQ.Ol  GO  TO  91 
00  83  JQ-itLYM 
S5=J0*Y2  . 

Ui=b2*CCOS(S9l+U*SIN(S9) ) 

00  84  JP  =  1 1  L X 
Jl=Jl*l 
P ( J 1  I =U1 
84  CONTINUE 
83  CONTINUE 
91  00  89  JQ=1 t  L  Y 
S5=  ( J0-. 81 *Y2 
Ul  =  Sl*(COS( S5I *U* S INI S  5 ) 1 
DO  86  JP  =  1  »  L  X  M 
J1=J1*1 
PI  J1)=U1 
86  CONTINUE 
89  CONTINUE 


RETURN 

FND 
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IV.  DESCRIPTION  OF  THE  PROGRAM  TO  PLOT  PATTERNS 


This  program  plots  patterns  of  the  transmission  cross  section  per 
square  wavelength  t/X'1'  read  from  direct  access  data  set  6. 


Punched  card  data  is  read  according  to 


READ (1,22)  NTH,  NP 
22  FORMAT (2013) 

READ  (1 ,22)  (I.P  (K)  ,  K  =  1,  NP) 
READ (1 , 10) (SCL(K) ,  K  =  1,  NP) 
10  FORMAT (6E11. 4) 

2 

The  patterns  of  t/X  are  read  from  direct 
REWIND  6 
J1  =  NTH*NP 

READ(6) (TA(I) ,  I  -  1,  Jl) 


access  data  set  6  according  to 


TA(J  +  (K-1)*NTH)  is  the  value  of  t/X2  at  angle  (J-l) *tt/(NTH-1)  radians 
on  the  Kth  pattern.  Here,  J  «  1,2, ...NTH  and  K  =  1,2. ..NP.  If  LP(K)  -  0, 
the  Kth  pattern  is  not  plotted.  If  LP  (K)  4  0,  the  Kth  pattern  is  multi¬ 
plied  by  SCL(K)  and  then  plotted  in  inches. 


Minimum  allocations  are  given  by 


DIMENSION  LP(NP) ,  SCL(NP),  SN(NTH),  CS(NTH), 

TA(NP*NTH) ,  X(NTH) ,  Y (NTH) 

DO  loop  15  plots  the  Kth  pattern  if  LP(K)  40.  DO  loop  16  puts 
tick  marks  on  the  vertical  axis  drawn  by  statement  25.  DO  loop  18  puts 
tick  marks  on  the  horizontal  axis  drawn  by  statement  26.  Statement  27 
plots  the  pattern  whose  horizontal  and  vertical  coordinates  have  been 
stored  in  X  and  Y  by  DO  loop  20. 
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C  LISTING  OF  THE  PROGRAM  TO  PLOT  PATTERNS 

C 

//  FXFC  FORTGCLG 
//FORT.SYSIN  00  * 

DIMENSION  LPI25),SCL(25),XX(4|,YY<4),SN(73),CS(73),TA<n68) 
DI MENS  ION  X(73I,Y(73) 

RF  AD  I  I  *  22 )  NTH, NP 

22  FOR  MAT ( 201 3  I 
WRITF(3,11)  NTH.NP 

11  FORMAT ( •  NTH  NP*/1X,2I3) 

RF AD( 1,22) ( LP(K) ,K=1,NP I 
WP I TE ( 3 , 23  I  I L  P ( K I  ,K=I,NP> 

23  FORMAT { •  LP*/I IX, 20! 3) ) 

RF ADII ,10M SCLIK) ,K=1,NP> 

10  FORMAT (6E1 I .4) 

WRITF (3,24) ISCLIK ) ,K=l,NP» 

24  FORMAT!*  SCL* /{  IX ,6Fil .41) 

XXIII  =  1. 

YY ( 1  )  *  1  • 

XXI 21=1. 

YY ( 2 ) =9. 

XX  I  3 ) =  1  . 

YY ( 3) =5. 

XX  (  4  )  =  5  . 

YY  (  4  )  =  5. 

PI =3. 141593 
DTH=P  I/INTH-l  ) 

00  19  JM.NTH 
ANG= ( J-l J*0TH 
SN ( J ) =  S ! NI ANG  > 

CS I J )=COS! ANG) 

19  CONTINUE 
CALL  PLOTID 
REWIND  6 

12  Jl=NTH*NP 
READI6MTAI  I),l=l,Jl) 

WR  I T  E  I  3 , 1 4  )  TAin 

14  FORMAT (•  TA=* , Ell .4) 

J1  =  0 

DO  15  K=l,NP 
IF(LPIK) .EO.O)  GO  TO  21 

25  CALL  L  INEI XXI 1 ) ,YY( 1 ) , 2,1 .0-0) 

S3  =  9  . 

DO  16  J=  1 , 9 

17  CALL  SYMBOL! 1., S3, .14, 13, 90., -I  I 
S3=  S  3-  l. 

16  CONTINUE 

26  CALL  L  INEI XX(3) ,YY<3) ,2 ,1 ,0,0) 

S 1  =  5  . 

DO  18  J= 1 , 4 

CALL  SYMBOL (SI, 5. ,.14, 13,0., -1) 

Sl*Sl-l. 

18  CONTINUE 
00  20  J= 1 , NTH 
J2  =  J 1 ♦ J 

Sl*T  A ( J2 ) *SCL I K ) 

X ( J ) = 1 .♦Sl*SN( J ) 

Y(  J)  =  5.*-Sl*CS(JI 

20  CONTINUE 

27  CALL  L INFIX ( 1) ,Y(  1 ) ,NTH,l ,0,0) 
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